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Abstract 

We investigate the integrated cross section for forward J/ip pro- 
duction in proton-proton and pion-proton collisions at centre of mass 
energies upto 60 GeV. We find that colour-octet contributions to the 
cross section are crucially important, and their inclusion produces a 
very good description of the data. Since the values of non-perturbative 
matrix elements are fixed by other experiments, these predictions are 
parameter- free. 



The production of bound states of heavy quarks has been theoretically 
studied in terms of the colour-singlet model [Q. In this approach, the 
dominant contribution to quarkonium production is expected to come from 
quark-antiquark or gluon-gluon fusion, leading to the formation of a heavy- 
quark pair in a colour- singlet state with the correct spin, parity and charge- 
conjugation assignments projected out. By using radial wave-functions or 
their derivatives derived from other measurements as inputs, the colour- 
singlet model makes definite predictions for the production cross sections 
of specific quarkonium resonances. 

While this model gives a reasonable description of large-p^ J/ ip produc- 
tion at ISR energies, it completely fails to explain the Pt integrated cross 
sections at fixed-target and ISR energies. This failure is easily traced back 
to the model assumptions. The requirement that a colour-singlet state with 
the correct quantum numbers be produced at the hard vertex imposes rather 
stringent constraints (coming from C-invariance) and rules out the zero-p-r 
production of the 5'-state quarkonia at lowest order. Thus, to leading order, 
the inclusive cross section for J/ip production is obtained only from direct 
production of Xc followed by their decay to J/iIj. 

Recent improvements in the theory of quarkonium production have been 
provoked by the anomalously large cross section for J/ ip production at large 
Pt measured recently in the CDF experiment at the Tevatron 0. For val- 
ues of Pt much larger than the charm quark mass, fragmentation of gluons 
and charm quarks become dominant The CDF data has been success- 
fully explained by taking into account both the fusion and fragmentation 
contributions [Q. 

The analyses of the CDF data have brought out another limitation of 
the colour-singlet approach, i.e., the relative velocity, f, between the heavy 
quarks in the bound state is ignored. However, v is often not negligible in 
quarkonium states and corrections of order v need to be taken into account. 
Starting from a non-relativistic QCD Lagrangian, a systematic analysis us- 
ing the factorisation method has been recently carried out 0]. In this for- 
mulation, the quarkonium wave-function admits of a systematic expansion 
in powers of v in terms of Fock-space components. For example, the wave- 
functions for the P-state charmonia have the conventional colour-singlet P- 
state component at leading order, but there exist additional contributions at 



1 



non-leading order in v, involving octet S'-state components; i.e., 

\Xj) = Oil)\QQ['Pi'^]) + Oiv)\QQ['S[%) + ... (1) 

While the octet components are necessary for a consistent perturbative de- 
scription of the P-states 0, for S'-states it is required in order to explain the 
anomalously large ip' direct J/ip production at the Tevatron [§]. 

As a consequence of factorisation, each matrix element is split into two 
parts. The short-distance parts for the colour-octet process are calculable, 
but the long-distance non-perturbative octet matrix elements (the analogues 
of the wavefunctions in the colour- singlet model) are not. They have been 
obtained by a fit to the Tevatron data More recently, the colour-octet 
processes which contribute to photoproduction of J/ip have been studied 
Using both the fixed-target and HERA data on photoproduction, a different 
linear combination of the same octet matrix-elements has been fitted. 

The time seems ripe for a reappraisal of the mass of data on hadropro- 
duction of J/ip in proton-nucleon and pion-nucleon collisions at centre of 
mass energies, v^, upto about 60 GeV (see [0] for a compilation of data). 
In this paper we carry out this investigation. We find that the low-energy 
data on forward {xp > 0) integrated hadroproduction cross sections are well 
described by the inclusion of octet components. These are parameter-free 
predictions, because the values of the non-perturbative matrix elements are 
derived from other experiments. 

In this model, a. J/ip can be obtained either directly, or through radiative 
decays of Xc states. Any of these can be produced through a colour-singlet 
channel or through the octet channel. Quantum number arguments show 
that the cross sections for direct J/ip or xi production in the singlet channel 
are zero (to lowest order in as)- As a result, we can write the cross section 



where BR^j denotes the branching ratio for the decay of Xj i^ito J/ip. The 
individual terms are all known. We write down the forward integrated cross 
sections below. 
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The singlet channel cross sections |1|] are given in terms of the strong 
coupling Os and the charm quark mass m by 

The quantity G is defined in terms of the gluon momentum densities in 
the projectile (P) and target (T) by 

G = [ —gp{x)gr{T/x). (4) 

X 

The lower limit of the integral, y/r = 2m/ ^/S, is written for the forward 
cross section. If all values of Xp were to be considered, the limit would have 
changed to r. 

The two matrix elements are related to the derivative of the wavefunction 
at the origin by 

orcpj)) = ^-^^\R'm'. (5) 



The octet cross section for J/^jj has been calculated recently [|Tl[]. It is 
given by 



2^3 



(6) 

The quantity Q is defined in terms of the quark momentum densities in the 
projectile and target by 

Q = T. -lii^Mir/x) + (P - T), (7) 

J- J \/r X 

where the sum is over fiavours. The combination of non-perturbative matrix 
elements in the coefficient of G is precisely that needed for photoproduction 

i- 

There are also the octet contributions for the Xc states. These can be eas- 
ily obtained by noting the following facts. The leading contribution comes 
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from the ^Si colour-octet states. The gluon fusion contribution to this chan- 
nel vanishes Hence, the only cross section required is 



The singlet matrix elements of eq. are extracted from data on hadronic 



decays of the Xc states. For a recent compilation see ||T^ . The colour-octet 
matrix elements {Og^{^Si)) and (^Og^^ {^Si)'^ have been extracted from the 
hadroproduction rates at the Tevatron. The different resonances are identi- 
fied in these experiments, allowing separate fits to these matrix elements. We 
use the tabulation of ||l3|. We require all the matrix elements (Cg"'(^S'i)). 



These are related by 

[ores,]) = i2J+l)(ori's^)). (9) 



The remaining combination of matrix elements has been extracted from pho- 
toproduction data 0. We use this value 

)J/i^(iQA\ _L JL _L ^ /nJ/-^(3TDA\\ _ o.020±0.001 GeV^ 

(10) 

in the calculations reported here. 

A different linear combination of the same matrix elements enters the 
large-pr hadroproduction cross sections relevant to the Tevatron experiments. 
A consistency check between these two linear combinations can be made 
if one assumes a relation between (O'l^'^ {^P2)^ and (ol^ Pq)^ . Such a 
relation reduces the number of unknown matrix elements by one and permits 
evaluation of the individual matrix elements. With the assumption that the 
first matrix element is 5 times the second, the Tevatron and photoproduction 
fits together yield a neg ative value of (C8/'^(3Po)) i- This may just indicate 
that the assumption has to be modified. The issue is, however, irrelevant to 
our computation, since we require precisely the linear combination in eq. 

We use the MRS D— ' and the GRV LO sets of parton densities for the 
proton. Since we are interested in values of \fS < 60 GeV, the lowest values 
of X entering the parton distributions are y/r ~ 0.05. For such large values of 
X the differences between structure function parametrisations for the proton 
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are rather small. The pion structure functions are not known as accurately. 
We have used the parametrisations SMRS 1, SMRS 3 and GRV. All the 
structure functions are taken from the PDFLIB package in . 

As always in lowest order QCD cross sections, there is a choice of scale to 
be considered. In addition, the charm quark mass, m is not precisely known. 



From open charm production, the limits are 1.2 GeV < m < 1.8 GeV 



Our results are shown for the two choices m = 1.7 GeV and a scale of 2m, 
as well as m = 1.6 GeV and a scale of 4m. 

With the inputs discussed here, we find that the y/S dependence of the 
integrated forward J/ijj production rates, for both pp and vrp collisions, are 
described rather well by the model (see Figure [^). We would like to emphasise 
that there are no free parameters in this calculation. 

We find that the total J/ip hadroproduction rate in pp collision is dom- 
inated by the gluon-gluon fusion part of the colour-octet direct J/ip cross 
section. Although xo is produced copiously, its contribution to the J/ip cross 
section is suppressed by a small branching ratio. The contribution of the xi 
is small because the singlet contribution is zero and the octet is small. The 
X2 state gives the major contribution among the P-channel resonances, but 
nevertheless gives less than 5% of the inclusive forward cross section. In up 
collisions the P-channel contributes slightly more to the total cross section. 
This is due to the enhancement of the qq initiated octet x contributions. 

Since the cross section is dominated by the gluon fusion process giving a 
colour-octet contribution to direct J/ip, the dominant uncertainty in the non- 
perturbative matrix elements comes from the linear combination in eq. ([T0|). 
This gives a 5% uncertainty in the predictions. Although the uncertainty in 
the non-perturbative matrix elements for colour-octet contributions to the x 
states are larger, they have a much smaller effect on the result. Similiarly, 
the assumption stated in eq. (^) has very little impact on the result. 

Larger uncertainties come from the choice of scale in QCD and the charm 
quark mass. The scale uncertainty can only be controlled by taking into 
account the order a| contributions to the cross section. When this is done, 
it might be possible to control the uncertainty in m. 

In conclusion, the world data on forward inclusive J/ip hadroproduction 
cross sections are well described by the inclusion of octet contributions. The 
calculation is completely parameter-free since the non-perturbative inputs 
are fixed from other experiments. Our results on other aspects of the phe- 
nomenology, such as rapidity and Xp distributions, will be reported elsewhere. 
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Figure 1: The colour-octet model predictions for integrated forward J/%Ij 
hadroproduction cross sections as a function of the CM energy. For pp colli- 
sions, the two curves for each choice of m and scale Q are for the structure 
functions MRS D— ' and GRV LO. For 7rp coUisions the three sets of structure 
functions are MRS D— ' for proton and SMRS 1 for tt, MRS D— ' for proton 
and SMRS 3 for tt, GRV LO for both. Note that the colour-singlet model 
predictions lie far below the data. 
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